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■ Abstract. We present a study of the shear modulus of the crystalline color superconducting phase 

^ ^ I of quark matter, showing that this phase of dense, but not asymptotically dense, quark matter 

Q . responds to shear stress as a very rigid solid. This phase is characterized by a gap parameter A that 

' is periodically modulated in space and therefore spontaneously breaks translational invariance. We 

derive the effective action for the phonon fields that describe space- and time-dependent fluctuations 
of the crystal structure formed by A, and obtain the shear modulus from the coefficients of the spatial 
derivative terms. Within a Ginzburg-Landau approximation, we find shear moduli which are 20 to 
' (— I ' ] 1000 times larger than those of neutron star crusts. This phase of matter is thus more rigid than any 

known material in the universe, but at the same time the crystalline color superconducting phase is 
also superfluid. These properties raise the possibility that the presence of this phase within neutron 
stars may have distinct implications for their phenomenology. For example, (some) pulsar glitches 
may originate in crystalline superconducting neutron star cores. 
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m : INTRODUCTION 

on ■ At the large baryon number densities found at the cores of neutron stars, quarks may not 

be confined into well-defined hadrons. If quarks are deconfined, since the interaction 
between quarks in the color antisymmetric channel is attractive, quarks near their Fermi 
I surfaces will tend to form Cooper pairs. Because the temperature of neutron stars is 

O ■ much smaller than the typical pairing energy, if quark matter exists within neutron stars, 

it must be in some color superconducting phase yj]. 
; At asymptotically large densities where the masses of the u, d and s quarks can all 

■ be neglected, quark matter exists in the CFL phase |0] in which quarks of all three 

colors and all three flavors form Cooper pairs with zero total momentum. The diquark 
condensate in the CFL phase is antisymmetric in color indices, driven by the attractive 
color interaction, and antisymmetric in the spin indices. Thus it must be antisymmetric 
in flavor indices, implying that quarks of a given flavor can only pair with quarks of the 
other two flavors. In the CFL phase, all fermionic excitations are gapped, with a gap 
parameter Aq ~ 10 - 100 MeV. 

However, at densities relevant for neutron star phenomenology, meaning quark chem- 
ical potentials at most /i ~ 500 MeV, the strange quark mass cannot be neglected. 
In neutral unpaired quark matter in weak equilibrium, Ms induces splitting between the 
Fermi surfaces for quarks of different flavor, which can be taken into account to low- 
est order in Mj/ji'^ by treating the quarks as if they were massless but with chemical 
potential splittings 5/i2 = (/ij, — /i.?)/^ and 5/i3 = (/i^ — /i;,)/2 given by 5/i2 = 5/i3 = 
SjJ. — M^/(8/i). Note that the splitting between unpaired Fermi surfaces increases with 



decreasing density. In the CFL phase, the Fermi momenta are not given by these optimal 
values for unpaired quark matter; instead, the system pays a free energy price oc 
to equalize all Fermi momenta and gains a pairing energy benefit Ag/i^. As a function 
of decreasing density, there comes a point (at which 8jl d^/A lH]) when the system 
can lower its energy by breaking pairs. Restricting the analysis to spatially homogeneous 
condensates, the phase that results when CFL Cooper pairs start to break is the gapless 
CFL (gCFL) phase [3]. However, this phase turns out to be "magnetically unstable" 
meaning that it is unstable to the formation of counter-propagating currents. If Aq is 
small enough that CFL pairing cannot survive all the way down to the /i at which quark 
matter is supplanted by nuclear matter, then the true ground state of intermediate density 
quark matter must have a lower free energy than that of the unstable gCFL phase. 

Crystalline color superconducting quark matter is a possible resolution of the mag- 
netic instability of the gCFL phase. Crystalline color superconductivity [l5|,@,|7|,[8l] is the 
QCD analogue of a form of non-BCS pairing first considered by Larkin, Ovchinnikov, 
Fulde and Ferrell [@] . This phase may be the ground state of matter in the intermedi- 
ate density regime in which quark matter is favored over nuclear matter but the Fermi 
surface separations oc / ji are large enough to disrupt CFL pairing. In a crystalline 
color superconducting phase, quarks whose Fermi surfaces are separated (as favored in 
the absence of pairing) nevertheless pair. Unlike in conventional BCS phases, the quarks 
in a pair do not have equal and opposite momenta (meaning that the Cooper pairs have 
net momentum) allowing both members in a pair to have momenta near their respective. 



separated, Fermi surfaces. Such phases do not suffer from the magnetic instability [|10|] . 
A particular crystalline phase is specified by sets of momentum vectors {q^}, meaning 
that Cooper pairs with total momentum 2q" form for each q? G {Q/}. All the q^'s have 
the same magnitude qi = \q"\ = r\djli, with r\ = 1.1997 ^ and 25/i/ the separation 
between the Fermi surfaces of quarks that pair. The directions of the vectors in {q^} 
must be determined to get structures with the smallest free energy. In position space, the 
condensate is 

( WmCr' Wjfi ) - E e/«^ ^liAi E exp (2/q^ ■ r) . (1) 
t qye{q/} 

This is antisymmetric in color (a,/3), spin, and flavor (/, j) (where (1, 2, 3) correspond 
to (m, d, s) respectively) indices and is thus a generalization of the CFL condensate to 
crystalline color superconductivity. For simplicity, Ai is set to 0, neglecting (ds) pairing 
because the d and s Fermi surfaces are twice as far apart from each other as each is from 
the intervening u Fermi surface. Hence, / can be taken to run over 2 and 3 only. {q2} 
and {qg} define the crystal structures of the {us) and {ud) condensates respectively. 

We will analyze crystalline color superconductivity in an NIL model, which gives in 
the mean field approximation an interaction term 

1 _ _^ 

^interaction = 2 "/^^(r) T +h.C., (2) 

where the proportionality constant in Eq. ([U) is conventionally chosen so that, 

A(r) = (Cr5)£e,„^e;,-,-A; £ exp(2/q^r). (3) 

1 q?e{q/} 



The authors of fsf] calculated the free energy Q. of several crystalline structures within 
the weak coupling (5/i,Ao <C /i) and Ginzburg-Landau (A <^ Aq) approximations, 
and found qualitative features that make a structure favorable. Two particular structures 
that possess these features, called CubeX and 2Cube45z and described below, have a 
lower Q. than any other crystal patterns that have been analyzed. One or the other or 
both of these two is favored over unpaired quark matter and the gapless CFL phase over 
the large range of densities given by [SJ 

2.9Ao < ^ < 10.4Ao . (4) 

For Ao = 25MeV and Ms = 250MeV, the window of densities over which one of these 
crystalline phases has the lowest free energy of any color superconducting phase extends 
from well below the ji at which nuclear matter supersedes quark matter to well above 
the highest /i ~ 500 MeV expected in the cores of neutron stars. The robustness of 
these phases is partly due to their having reasonably large gap parameters A, so large 
that the Ginzburg-Landau expansion parameter (A/5/i)^ can be around a tenth to a 
fourth meaning that this approximation is at the edge of its validity. Nevertheless, their 
impressive robustness over a large range of /i relevant for cores of neutron stars make 
it worth considering the phenomenological implications of their presence. We also note 
that in the case of a simpler crystal pattern, for which results have been obtained without 
making the Ginzburg-Landau approximation, this approximation is conservative in that 
it always underestimates both A and the condensation energy [7]. 

In the CubeX crystal structure, {qj} and {q3} each contain four unit vectors, 
with {q2} = {(l/y3)(±y2,0,±l)} and {qa} = {(lv^)(0,±y2,±l)}. In the 
2Cube45z crystal structure {qj} and {q^} each contain eight unit vectors, with {q2} = 

{(l/v^)(±l,±l,±l)} and {%} = {(1/V3)(±v^,0,±l)} U{(l/V3)(0,±v^,±l)}. 
For these structures, {q2} can be exchanged with {qa} by rigid rotations, ensuring that 
there are electrically neutral solutions of the gap equation with A2 = A3 = A [js'], a fact 
that we will use in the next Section. 



PHONONS 

The crystalline phases of color superconducting quark matter that we have described in 
the previous Section are unique among all forms of dense matter that may arise within 



neutron star cores in one respect: they are rigid UllH . They are not solids in the usual 
sense: the quarks are not fixed in place at the vertices of some crystal structure. Instead, 
these phases are in fact superfluid since the condensates all spontaneously break the 
f/(l)g symmetry corresponding to quark number. The diquark condensate, although 
spatially inhomogeneous, can carry supercurrents [[sl, [H]. And yet, we shall see that 



crystalline color superconductors are rigid solids with large shear moduli. It is the spatial 
modulation of the gap parameter that breaks translation invariance, and it is this pattern 
of modulation that is rigid. This novel form of rigidity may sound tenuous upon first 
hearing, but we shall present the effective Lagrangian that describes the phonons in 
the CubeX and 2Cube45z crystalline phases, whose lowest order coefficients have been 



calculated in the NJL model that we are employing fll']. We shall then extract the shear 
moduli from the phonon effective action, quantifying the rigidity and indicating the 
presence of transverse phonons. 

The shear moduli of a crystal may be extracted from the effective Lagrangian that 
describes phonons in the crystal, namely space- and time-varying displacements of the 
crystalline pattern UM- In the present context, we introduce displacement fields for the 
{ud), {us) and {ds) condensates by making the replacement 



q?e{q,} 



,2<.(r-u/(r)) 



(5) 



in ([3]). One way to obtain the effective action describing the dynamics of the displace- 
ment fields U/(r), including both its form and the values of its coefficients within the 
NJL model that we are employing, is to take the mean field NJL interaction given by 
Eq. (l2l), but with ([5]), and integrate out the fermion fields. Since the gapless fermions do 
not contribute to the shear modulus, one can integrate out the fermions completely for 
this calculation. As an aside, we mention that this is not true for the calculation of ther- 
mal or transport properties, where the gapless fermions do contribute and indeed may 
dominate, as is the case for the heat capacity [|l3] an d neutrino emissivity lll4ll . 
After integrating out the fermions, we obtain jlll] 



S[u] 



d^x ^ Ki 



X 



.qye{q/} 



vq?e{q,} 



(6) 



where m, n, v and w are spatial indices running over x, y and z and where we have 

defined Ky = ^^ij^p^. For Ai =0, A2 = A3 = A, and r\ ~ 1.1997, K2 = k^ = kc^ 

0.664jU^|A^|. 5[u] is the low energy effective action for phonons in any crystalline 
color superconducting phase, valid to second order in derivatives, to second order in 
the gap parameters A/ and to second order in the phonon fields u/. Because we are 
interested in long wavelength, small amplitude, phonon excitations, expanding to second 
order in derivatives and in the phonon fields is satisfactory. The Ginzburg-Landau 
approximation is an expansion in (A/5/i)^, and so is not under quantitative control 
for the most favorable phases, as we have discussed. But, the main requirement from 
glitch phenomenology for the shear modulus is that it should be large, and given that we 
find much larger values than those obtained for conventional neutron star crusts, there 
is no great motivation to go to higher order to improve the precision. (A higher order 
calculation would include couplings between the different u/, meaning they could no 
longer be treated independently.) 



According to the theory of elastic media IllSh . the shear moduli can be extracted from 
the phonon effective action. Introducing the strain tensor 



1 /(9m™ du] 



(7) 



we then wish to compare the action ^ to 



5[u] = Ud'x f i:i:pr('?o<)(<?o<) - LLAr""'^^ , (s) 

^ J \ r m I nin I 



I m I 

vw 



which is the general form of the action quadratic in displacement fields and which 
defines the elastic modulus tensor X^^'"''^ . In this case, the stress tensor (in general the 
derivative of the potential energy with respect to s'f^') is given by 

Oi =Ai Sj . (9) 

The diagonal components of o are proportional to the compression exerted on the system 
and are therefore related to the bulk modulus of the crystalline color superconducting 
quark matter. Since unpaired quark matter has a pressure ~ /i^, it gives a contribution 
to the bulk modulus that completely overwhelms the contribution from the condensation 
into a crystalline phase, which is of order /i^A^. We shall therefore not calculate the 
bulk modulus. On the other hand, the response to shear stress arises only because of the 
presence of the crystalline condensate. The shear modulus is defined as follows. Imagine 
exerting a static external stress (J/ having only an off-diagonal component, meaning 
^mv ^ for a pair of space directions m ^ v, and all the other components of o are zero. 
The system will respond with a strain s"^. The shear modulus in the mv plane is then 

n-""' 1 

,,mv / T mvmv /■! m 

' =2^ = 2^ ' 

where the indices m and v are not summed. For a general quadratic potential with o^^ 
given by Q, v""' simplifies partially but the full simplification given by the last equality 
in (flOl) only arises for special cases in which the only nonzero entries in with 
m^v are the A"'^™^ entries, as is the case for all the crystal structures that we consider. 

For a given crystal structure, upon evaluating the sums in Q and then using the 
definition (|7]) to compare © to ([8]), we can extract expressions for the A tensor and 



thence for the shear moduli. This analysis, described in detail in llllll . shows that in the 
CubeX phase 





V2=yK|^ I , V3 = — K| 1 I , (11) 
while in the 2Cube45z phase 

V2 = ^k| 1 1 I , V3 = ^k( 1 I . (12) 





We shall see in the next Section that it is relevant to check that both these crystals have 
enough nonzero entries in their shear moduli V/ that if there are rotational vortices are 



pinned within them, a force seeking to move such a vortex is opposed by the rigidity of 
the crystal structure described by one or more of the nonzero entries in the V/. This is 



demonstrated in Ullll . 

We see that all the nonzero shear moduli of both the CubeX and 2Cube45z crystalline 
color superconducting phases turn out to take on the same value, 



where /i is expected to lie between 350 and SOOMeV and A may be taken to lie between 
5 and 25MeV to obtain numerical estimates. 

From (fT3l) we first of all see that the shear modulus is in no way suppressed relative 
to the scale ji^A^ that could have been guessed on dimensional grounds. And, second, 
we discover that a quark matter core in a crystalline color superconducting phase is 20 
to 1000 times more rigid than the crust of a conventional neutron star [[l6ll . Finally, one 
can extract the phonon dispersion relations from the effective action The transverse 
phonons, whose restoring force is provided by the shear modulus turn out to have 
direction-dependent velocities that are typically a substantial fraction of the speed of 
light, being given by and a/2/3 UM- This is another way of seeing that this 

superfluid phase of matter is rigid indeed. 



RIGID QUARK MATTER AND PULSAR GLITCHES 

The existence of a rigid crystalline color superconducting core within neutron stars may 
have a variety of observable consequences. For example, if some agency (like magnetic 
fields not aligned with the rotation axis) could maintain the rigid core in a shape 
that has a nonzero quadrupole moment, gravity waves would be emitted. The LIGO 
non-detection of such gravity waves from nearby neutron stars ITtIi already limits the 
possibility that they have rigid cores that are deformed to the maximum extent allowed 
by the shear modulus (fT3]) . upon assuming a range of breaking strains, and this constraint 



will tighten as LIGO continues to run [11811 . Perhaps the most exciting implication of a 
rigid core, however, is the possibility that (some) pulsar "glitches" could originate deep 
within a neutron star, in its quark matter core. 

A spinning neutron star observed as a pulsar gradually spins down as it loses rotational 
energy to electromagnetic radiation. But, every once in a while the angular velocity 
at the crust of the star is observed to increase suddenly in a dramatic event called a 
glitch. The standard explanation lfl9ll requires the presence of a superfluid in some 
region of the star which also features a rigid structure that can pin the vortices in the 
rotating superfluid and that does not easily deform when the vortices pinned to it are 
under tension. As a spinning pulsar slowly loses angular momentum over years, since 
the angular momentum of any superfluid component of the star is proportional to the 
density of vortices, the vortices "want" to move apart. However, if some vortices are 
pinned to a rigid structure and so do not move, after a time this superfluid component 
of the star is spinning faster than the rest of the star. When the "tension" built up in the 
array of pinned vortices reaches a critical value, there is a sudden "avalanche" in which 



vortices unpin, move outwards reducing the angular momentum of the superfluid, and 
then re-pin. As this superfluid suddenly loses angular momentum, the rest of the star, 
including in particular the surface whose angular velocity is observed, speeds up — a 
glitch. In the standard explanation of pulsar glitches, the necessary conditions are met 
in the inner crust of a neutron star which features a neutron superfluid coexisting with 
a rigid array of positively charged nuclei that may serve as vortex pinning sites. In very 
recent work. Link has questioned whether this scenario is viable because once neutron 
vortices are moving through the inner crust, as must happen during a glitch, they are so 
resistant to bending that they may never re-pin [lioll . Link concludes that we do not have 
an understanding of any dynamics that could lead to the re-pinning of moving vortices 
in the crust, and hence that we do not currently understand the origin of glitches as a 
crustal phenomenon. 

By virtue of being simultaneously superfluids and rigid solids, the crystalline phases 
of quark matter provide all the necessary conditions to be the locus in which (some) 
pulsar glitches originate. Their shear moduli (fT3l) make them more than rigid enough for 
glitches to originate within them. The crystalline phases are at the same time superfluid, 
and it is reasonable to expect that the superfluid vortices that result when a neutron 
star with such a core rotates have lower free energy if they are centered along the 
intersections of the nodal planes of the underlying crystal structure, i.e. along lines 
along which the condensate already vanishes in the absence of a rotational vortex. A 
crude estimate of the pinning force on vortices within crystalline color superconducting 
quark matter indicates that it is sufficient [|lll] . So, the basic requirements for superfluid 
vortices pinning to a rigid structure are all present. The central questions that remain 
to be addressed are the explicit construction of vortices in the crystalline phase and 
the calculation of their pinning force, as well as the calculation of the timescale over 
which sudden changes in the angular momentum of the core are communicated to the 
(observed) surface, presumably via the common electron fluid or via magnetic stresses. 

Much theoretical work remains before the hypothesis that pulsar glitches originate 
within a crystalline color superconducting neutron star core is developed fully enough 
to allow it to confront data on the magnitudes, relaxation timescales, and repeat rates 
that characterize glitches. Nevertheless, this hypothesis offers one immediate advantage 
over the conventional scenario that relied on vortex pinning in the neutron star crust. It 
is impossible for a neutron star anywhere within which rotational vortices are pinned to 
precess on ~ year time scales 12111. and yet there is now evidence that several pulsars 
are precessing [|22l] . Since all neutron stars have crusts, the precession of any pulsar is 
inconsistent with the pinning of vortices within the crust, a requirement in the standard 
explanation of glitches. On the other hand, perhaps not all neutron stars have crystalline 
quark matter cores — for example, perhaps the lightest neutron stars have nuclear matter 
cores. Then, if vortices are never pinned in the crust but are pinned within a crystalline 
quark matter core, those neutron stars that do have a crystalline quark matter core can 
glitch but cannot precess while those that don't can precess but cannot glitch. 
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